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Abstract 

We report the performances of a 0.51 kg CdW04 scintillating bolometer to be used for future Double Beta Decay Experiments. 
The simultaneous read-out of the heat and the scintillation light allows to discriminate between different interacting particles aiming 
at the disentanglement and the reduction of background contribution, key issue for next generation experiments. We will describe 
the observed anticorrelation between the heat and the light signal and we will show how this feature can be used in order to increase 
the energy resolution of the bolometer over the entire energy spectrum, improving up to a factor 2.6 on the 2615 keV line of ^"**T1. 
The detector was tested in a 433 h background measurement that permitted to estimate extremely low internal trace contaminations 
of ^-^^Th and ^^^U. The light yield of yjfi, a's and neutrons is presented. Furthermore we developed a method in order to correctly 
evaluate the absolute thermal quenching factor of a-particles in scintillating bolometers. 
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1. Introduction 

Double Beta Decay (DBD) searches became of critical im- 
portance after the discovery of the neutrino oscillations and 
plenty of experiments are now in the construction phase and 
many others are in R&D phase |[T]|2]|3]|4l. The main challenges 
for all the different experimental techniques are the same [5 1: i) 
increase the active mass, ii) decrease the background, and iii) 
increase the energy resolution. 

Thermal bolometers are ideal detectors for this survey: crys- 
tals can be grown with different interesting DBD-emitters and, 
fundamental for next generation experiments, they show an ex- 
cellent energy resolution. 

The CUORICINO Experiment ||6l, constituted by an array of 
62 Te02 crystal bolometers, demonstrated not only the power 
of this technique but also that the main source of background 
for these detectors arises from surface contaminations of ra- 
dioactive ff-emitters. Moreover simulations show that this con- 
tribution will largely dominate the expected background of the 
CUORE Experiment QEl in the region of interest, since there 
is no possibility to separate this background from the two DBD 
electrons. The natural way to discriminate this background is to 
use a scintillating bolometer [9|. In such a device the simulta- 
neous and independent read out of the heat and the scintillation 
light permits to discriminate events due to y/yS, a and neutrons 
thanks to their different scintillation yield. Moreover if the crys- 
tal is based on a DBD emitter whose transition energy exceeds 
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the 2615 keV y-line of ^o^xi then the environmental background 
due to natural y's will decrease abruptly. 

CdW04 is an ideal candidate for such kind of detector: 

• it is a well established scintillator 

• '""Cd (7.5 % i.a) has a DBD transition at 2805 keV 

• the light yield (LY) is rather large 

• the radiopurity of this compound is "naturally" high 

Due to these favourable features this crystal compound was 
already used to perform a DBD experiment IfTOl using "stan- 
dard" Photomultipliers. A large mass experiment, based on 
enriched CdW04 crystals, readout by Photomultipliers, was 
also proposed ifTTIl . But this technique, limited by the modest 
achievable energy resolution, key point of future experiments, 
is no more pursued. 



2. Experimental details 

The dimensions of the CdW04 crystal tested are 4 cm diam- 
eter, 5 cm height. All the surfaces of the crystal are polished at 
optical grade. The crystal was tested as a standard scintillator at 
room temperature. It was wrapped with a reflecting sheet (3M 
Radiant Mirror film VM2000) and coupled with optical grease 
to a Photomultiplier (Hamamatzu R6233). The energy resolu- 
tion evaluated on the '"Cs line is 12.3 % FWHM. 
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Figure 1 : Setup of the detector. The two Cu columns, holding the two frames, 
are not visible from the chosen perspective. 



Our detector setup is schematized in Fig. [T] It is held by 
means of four L-shaped Teflon pieces fixed to the two cylindri- 
cal Cu frames. The frames are held together through two Cu 
columns. The crystal is surrounded (without being in thermal 
contact) by a 43 mm diameter reflecting sheet (3M VM2002). 
The Light Detector (LD) |T2] is constituted by a 36 mm diam- 
eter, 1 mm thick pure Ge crystal absorber. The surface of the 
Ge facing the crystal is further "darkened" through the depo- 
sition of a 600 A layer of SiOa in order to increase the light 
abortion. Furthermore since the diameter of the LD is smaller 
with respect to the one of the reflecting cavity, a "reflecting 
ring" has been mounted in order to decrease light losses. On 
the opposite face of the crystal a reflecting sheet is mounted. 
The temperature sensor of the CdW04 crystal is a 3x3x1 mm^ 
neutron transmutation doped Germanium, the same used in the 
CUORICINO experiment. The temperature sensor of the LD 
has smaller volume (3x1.5x0.4 mm^) in order to decrease its 
heat capacity, increasing therefore its thermal signal. A resis- 
tor of ~300 kQ, realized with a heavily doped meander on a 
3.5 mm^ silicon chip, is attached to each absorber and acts as 
a heater to stabilize the gain of the bolometer llT3l [T4l . The 
detectors were operated deep underground in the Gran Sasso 
National Laboratories in the CUORE R&D test cryostat. The 
details of the electronics and the cryogenic facility can be found 
elsewhere ifTSi [T6l flTl . 

The heat and light pulses, produced by a particle interacting 
in the CdW04 crystal and transduced in a voltage pulse by the 
NTD thermistors, are amplified and fed into a 16 bit NI 6225 
USB ADC unit. The entire waveform of each triggered voltage 
pulse is sampled and acquired. The time window has a width 
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Figure 2: Scatter plot Light vs. Heat obtained in a 96 h calibration using an 
external ^^^^Th source. In the inset the highlight of the ^os-pj [jj^g jj^ jjjg circle 
two events due to the internal a-decay of '^"W. The energy spectrum below 400 
keV is completely dominated by the "'Cd/3-decay, with a rate close to 0.4 Hz. 



of 256 ms sampled with 512 points. The trigger of the CdW04 
is software generated while the LD is automatically acquired in 
coincidence with the former. The amplitude and the shape of 
the voltage pulse is then determined by the off line analysis that 
makes use of the Optimal Filter technique. The energy calibra- 
tion of the CdW04 crystal is performed using y sources placed 
outside the cryostat. The Heat axis is calibrated attributing to 
each identified y peak the nominal energy of the line, as if all 
the energy is converted into heat. Consequently this calibration 
does not provide an absolute evaluation of the heat deposited in 
the crystal. The dependency of amplitude from energy is pa- 
rameterized with a second order polynomial in log(V) where V 
is the heat pulse amplitude. The three coefficients of the poly- 
nomial are fitted on calibration data (the heat pulse amplitude 
corresponding to the known y lines visible in the spectrum). 
The choice of such a function was established by means of sim- 
ulation studies based on a thermal model of the detectors. 

The energy calibration of the LD is obtained thanks to a weak 
^^Fe source placed close to the Ge that illuminates homoge- 
neously the face opposed to the CdW04 crystal. During the 
LD calibration its trigger is set independent from the one of the 
CdW04. The LD is calibrated using a simple linear function. 
The FWHM energy resolution of the LD, evaluated on the X 
doublet at 5.90 and 6.49 keV, is 480 eV FWHM. In the range 
of Oh-50 keV, i.e. the energy interval of the light signals, the 
response of the LD can be definitely assumed to be linear. 

Three sets of data have been collected with this device: two 
calibrations using ^^^Th and ^"K sources, a long background 
measurement (433 h) and a neutron measurement (8 h) done 
exposing the detector to an Am-Be source. 

3. Light-Heat scatter plot 

In this field the usual way to present the results is to draw 
the Light vs. Heat scatter plot. Here each event is identified 



by a point with abscissa equal to the heat signal (recorded by 
the CdW04 bolometer), and ordinate equal to the light signal 
(contemporary recorded by the LD). In the scatter plot, yip, a 
and neutrons give rise to separate bands, in virtue of their char- 
acteristic Light to Heat ratio (see Figs. |5] and |6]). This feature 
is the result of the different LY's characterizing these particles: 
yIP events belong to a distribution (see Fig.[3]l characterized by 
the 9py angle, while the a events are characterized by 9a. 

The scintillation Quenching Factor (QF) is defined as the ra- 
tio of the scintillating yield of an interacting particle {a, neu- 
tron, nucleus) with respect to the LY of a ylfi event at the same 
energy. From Fig. [3] we have QFQ.=tan(0„)/tan(00y). But this 
holds only in first approximation, as will be exposed in Sec. l5] 
and Sec. |6] 

Within each band, monochromatic events (i.e. those corre- 
sponding to the complete absorption of a monochromatic parti- 
cle in the CdW04 crystal) appear as sections of lines with neg- 
ative slope, showing a strong anticorrelation between Heat and 
Light. This is evident in the Light vs. Heat scatter plot ob- 
tained with a source of ^^^Th (Fig. pb where different y-lines 
are clearly visible. Despite the much lower statistics, the same 
feature is perceptible also for the a lines that appear in the back- 
ground spectrum (Fig.lTll. 

A simple model accounts for the observed pattern. The 
energy E of a monochromatic particle, fully absorbed in the 
CdW04 crystal, is divided into two channels: a fraction (L) is 
spent in the production of Light (photons) and a fraction (//) 
is spent in the production of Heat (phonons). In absence of 
"blind" channels (i.e. channels in which the energy deposited 
into the crystal is stored in some system that do not takes part 
in signal formation) the energy conservation requires: 



E = Enea, + Eugh, ^(l-k)E + kE^H + L 



(1) 



being k the fraction of the total energy that escapes the crystal in 
form of light. The last equation states not only the energy con- 
servation, but also a subtle distinction between the Heat (what 
is measured) and the energy definition. 

Monochromatic events should produce a spot in the scatter 
plot, with a size determined by the intrinsic energy resolution 
of the heat and light detectors. The observed spread along a 
negative slope line, however, can be accounted for assuming 
the existence of fluctuations in the H and L signal amplitudes. 
Since the energy has to be conserved the amplitude of the fluc- 
tuations (6L and SH) must compensate each other: SL--6li. In 
other words the fluctuations are anti-correlated between each 
other and produce the observed pattern. They are visible be- 
cause their amplitude is much larger than the intrinsic resolu- 
tion of the detectors. 

We can devise two mechanisms that introduce a fluctuation 
in the L/H ratio: i) the statistical Poissonian fluctuation in the 
emitted light: the more the energy spent in light production, the 
less the energy spent in heat; ii) any variation of the LY as those 
due to position dependent effects (inhomogeneities and defects 
of the crystals that modify the light emission or self absorption 
of scintillation photons). Regardless the mechanism, we will 
demonstrate, in the next section, how this anticorrelation can 
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Figure 3: Light vs. Heat scheme for the interpretation of the energy correction 
method. The blue (dark) spots represent the y//3 monocliromatic events, while 
the red (light) represents an a. 



be used in order to increase (correct) the energy resolution of 
the detector over the entire energy spectrum. 

Fig. [3] tries to summarize the model here discussed. The 
dark spots represent monochromatic energy deposition in the 
CdW04 crystal, as those observed during a y calibration. For 
each of them a negative slope (iso-energy) line is tracked. The 
intercept on the Heat axis (H^'^'""") is the Heat that would cor- 
respond to a full heat conversion of the deposited energy (in 
absence of light emission or in the case in which all the light is 
absorbed by the crystal itself). In this simple model y/j6's and 
a's releasing the same energy within the CdW04 crystal lie on 
the same iso-energy line. But due to the larger LY, y//3 events 
will convert less heat within the crystal with respect to an a 
particle releasing the same energy: part of the energy escapes 
the crystal in form of photons. This feature can be easily de- 
duced from Fig.l3] the ai particle release the same total energy 
of the 72 particle (they belong to the same iso-energy line) but 
they show a different position on the Heat axis. In fact the o-i 
shows the same heat of the yi particle that releases a larger en- 
ergy within (but not "into") the crystal. The "usual" picture (in 
which y/p and a of the same energy show the same heat signal) 
corresponds to the limit in which the light emitted is negligi- 
ble (or completely re-absorbed by the crystal): in this case the 
distinction between H and E becomes meaningless. 

It is important to remark that up to now it was assumed to 
be able to measure the absolute values of energy converted into 
heat and light. 

Actually, the measured experimental heat and light signals 
can be written as 



i/ = a(l -k)E 
L = pekE 



(2) 



being a and p the absolute calibration factors for the Heat and 
Light axis. The factor e, instead, takes into account the overall 
light collection efficiency. 



Under the condition that a and /3 are not depending on the 
energy, and considering that the conservation of energy requires 
that 6[kE]^- 6[(l - k)El Eq. |2]impHes that 
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-tan(9fi,) 



(3) 



The last equation states that the effect of miss-caUbrations of 
the Heat and Light axis simply implies a variation of the slope 
of monochromatic lines, which is exactly what we observe ex- 
perimentally. Also the light collection efficiency plays a funda- 
mental role in the evaluation of the slope, as discussed in the 
next section. 

It is important to note that, unlike the Heat channel, the LD 
energy calibration is more delicate. The energy calibration on 
this detector is performed using ionizing X-Rays from the ^^Fe 
source. A priori the thermal signal that arises by the absorption 
of 1 keV of photons at the boundaries of the Ge crystal (more- 
over covered by the SiOi layer) could give a different thermal 
signal with respect to a 1 keV ionizing energy. In any case this 
effect still preserves proportionality, so that it can be considered 
within the factor /3. 
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Figure 4: FWHM energy resolutions evaluated for different 6 angles. For each 
angle a spectrum is produced and the corresponding peaks are fitted. The curves 
through the data points are only to guide eyes. For two peaks the corresponding 
error on the FWHM is not plotted for better clarity. 



4. Light-Heat anticorrelation: energy resolution 

The energy anticorrelation between the two signals of 
a double readout system was already observed in ioniza- 
tion/scintillation detectors 1 1 8 1 as well as in heat/scintillation 
bolometers |,19j with a device very similar to the one presented 
here. In the former case this anticorrelation was demonstrated 
to improve the energy resolution by a factor close to ~ 25 % 
(evaluated @570 keV). In paper 1 19 1 two different y-lines were 
studied at rather low energy. An evident anticorrelation was 
found but, however, the correlation factor of the two lines was 
found to strongly depend on the energy. 

The energy correction procedure is graphically explained in 
Fig. [3] Here blue (dark) and red (light) spots represent y/jS 
and a produced by monochromatic events. In order to obtain a 
spectrum with improved energy resolution we have to combine 
the heat and light values of each single event in an appropriate 
way. This can be obtained in a natural way performing a simple 
SU(1) rotation in the scatter plane an projecting the points on 
the Heatfl axis: 



Heatg = Heat cosfl -i- Light sinO 
Lighe = -Heat sin6' + Light cosO 



(4) 



The value of can be evaluated in two different ways: i) it 
can be "optimized" in order to obtain the best energy resolu- 
tion on the Heatfl axis (energy minimization); ii) by fitting the 
monochromatic spots in the scatter plot with a negative slope 
line, evaluating the 0/,-, angle for each single distribution (for 
construction we have 6 = n/l - Ofi,) 

It is clear from this scheme that some assumptions should be 
satisfied in order for this procedure to work properly: 

• the slope of each monochromatic spot (i.e. the 9fi, angle) 
or the 9 angle minimizing the FWHM has to be the same 
within all the lines; 



• the LY of the different class of interacting particles (a's 
and y/yS) has to be independent from energy. 

The first condition ensures that the rotation minimizes the 
FWHM on the whole energy spectrum. It has to be noted that if 
the first condition is not verified, this technique still provides a 
useful tool for resolution optimization. In case of different 0/„ 
values the rotation angle will be chosen in order to reach the 
best performances in the region of more interest for the physics 
(e.g. the region where the DBD peak should appear). 

The second condition ensures that the projected spectrum 
keeps linearity and energy calibration. This point is less triv- 
ial and will be treated in more detail in Sec.|5] 

We applied both the methods in order to evaluate the rota- 
tion angle. Using the energy minimization method, the value of 
the rotation 9 angle is allowed to vary over a wide range (0-80 
degrees) and for each angle a spectrum is produced projecting 
the data on the Heatg axis. Finally the FWHM energy resolu- 
tion of the more intense peaks is evaluated by a fit (we use an 
asymmetric Gaussian reporting then the average FWHM). 

We applied this method not only to the cahbration spectrum 
but also to the background measurement. In this way we were 
able to study also the the 88 keV y-line of '"''Cd (accidentally 
contamination of our CdW04 sample) and the a peaks due to 
internal trace contaminations: in this latter case we observe 
a line at the full energy of the decay (Q=a-i-nuclear recoil). 
The internal a-lines observed are 2^'^U-a(Q=4858 keV), '^''W-a 
(Q=2516 keV) and ^^^\J-a (Q=4270 keV). 

In Fig. p] we report the FWHM energy resolutions evaluated 
on the most intense lines obtained varying the rotation angle 9 
and projecting on Heate . The energy resolutions of the peaks 
are evaluated for each single rotated spectrum. The minimum 
of the energy resolution occurs within 9-{67.5 + 0.5) deg for all 
the lines, except for the 88 keV Electron Capture internal line 



of lo^Cd. 

In fact the minimum of the energy resolution for this Hne, 
1.03 ± 0.06 keV FWHM is obtained at 0= 60 deg while at 67.5 
deg it becomes 1.11 ±0.06 keV. These two values seem rather 
similar, but from Fig. [4] it turns out that the resolution curve 
for this peak shows a different behaviour with respect to all the 
other lines. We do not have an explanation for this enhancement 
in the anticorrelation at low energies, but the same effect was 
akeady observed in ||19|. 

Each monochromatic events distribution in the scatter plot 
was also fitted with a linear function, as previously discussed. 

The 9fi, angles as well as the most relevant parameters of this 
analysis are summarized in Tab.fTl The values of the 0/,, angles 
for the Q-'s and y/jS's lines reported here are compatible within 
the experimental errors within the two groups (with the excep- 
tion of the 88 keV line). However, with respect to the y/jS's, 
the angular coefficient calculated by fitting linearly the single 
lines gives a result that is systematically slightly larger with re- 
spect to the value (67.5+0.5) obtained through the energy min- 
imization. This systematic slightly difference could arise from 
the fact that while the bidimensional fit is model independent, 
the energy minimization depends on the fitting procedure of the 
peaks. Moreover it should be pointed out that the quoted error 
on 6fi, don't include systematic effects (introduced, for exam- 
ple, by the continuum background) whose values could domi- 
nate the total error 

The values of the 0/„ angles reported in Tab.fllseems to show 
a slightly difference between a and y//3. However it has to 
be pointed out that the quoted error on the 9fi, for the a par- 
ticles has a (further) systematic error that cannot be easily eval- 
uated. The point is that the a events were collected during the 
long background measurement. During this period the CdW04 
crystal was slowly cooling down (the LD, on the contrary, was 
rather stable during the same period). This drift, as explained 
in Sec. l2] is usually corrected through the Heater pulse. If the 
temperature drift is rather small (< 50 //K) then the correction 
is independent from the energy. On the contrary, if the drift 
is larger (in this case the drift was « 0.3 mK ) then the shape 
of the signal changes appreciably. In such cases the correction 
starts to depend on the energy. As a consequence a small drift 
in function of the time (or the temperature) becomes apprecia- 
ble. Moreover this drift depends on the energy. So, in order to 
evaluate the second, third and last column in Tab.[T]each a peak 
was "self-stabilized". This can introduce an uncontrolled sys- 
tematics since we are somehow forcing a distribution to obtain 
a minimum. On the other hand, a check was made on the weak 
2615 keV y-line present in the background spectrum. As in the 
case of the c's, we had to self-correct the drift. The obtained 
Ofi, value is 71.5±0.4, rather consistent with the a background. 
For the following we will assume Qf!^ — ffj-^. The ^■'^Th and 
^°K calibrations, as well as the 88 keV analysis, on the other 
hand, were made at the end of the measurement in a rather sta- 
ble working temperature. 

As a final remark it has to be noticed that the projection of 
the events on the Heatg axis introduces an Heat (energy) nor- 
malization factor K that will be different for different classes of 



particles, namely Ky^„=cos0(l-i-tan6' tanff^^Q,). This feature is 
rather evident in Fig. [3] in the Heat axis we have H(al)=H(yl) 
while in the rotated Heato axis we have H(a'l)=H(y2). This has 
two consequences: the first, obvious, is that the intercalibration 
between a and y/jS is different within the Heat spectrum and 
the Heatg spectrum. The second is that the same consideration 
holds within the same class of events: if the LY of a particles 
(for example) depends slightly on energy, then the Unearity of 
the projected spectrum will change. 

5. Light Yield and Scintillation Quenching Factor 

In Tab. [T] we also present the values of the LY, determined 
on the different monochromatic lines. We define the LY as the 
energy released in the LD (in keV) for a nominal energy depo- 
sition of 1 MeV in the scintillating crystal. 



Table 1 : Table with the main parameters of the detector. The theoretical FWHM 
energy resolution of the CdW04 crystal, evaluated through Optimal Filtering, 
is 0.8 keV. The LY of the a-particles (last three lines) is evaluated for the overall 
Q-value of the decay. The coirect value (assuming negligible scintillation for 
the nucleus recoil) is x: 2% larger. tSee text. 



Line 


FWHMe=o 


FWHMe=67.5 


LY 


n/l - Ofi, 


[keV] 


[keV] 


[keV] 


[keV/MeV] 


[deg] 


88 


1.4+0.08 


1.11+0.06 


17.3+0.2 


61.6+0.9 


583 


5.2+0.4 


2.1±0.2 


17.0+0.35 


70.8+0.4 


911 


8.4±0.5 


3.05+0.2 


17.5+0.2 


70.6+0.3 


1461 


11.7+0.5 


3.27+0.2 


17.45+0.15 


70.4+0.3 


2615 


16.5+0.5 


6.25+0.22 


17.6+0.1 


69.8+0.2 


2516 


5.2+1.6 t 


4.3+1.3 t 


2.74+0.08 


73.0+0.2 t 


4270 


7.0±1.5t 


4.4+1.0 t 


3.29+0.07 


72.3+0.2 t 


4858 


7.8±1.2t 


4.8+0.7 t 


3.44+0.08 


72.7+0.2 t 



The table shows that the LY of the y/jS events is, within the 
error, the same for all the observed lines (including the 88 keV 
line that shows a different resolution curve in Fig. |4|. On the 
other hand a small, but rather clear, energy dependence is evi- 
dent for a particles. This effect can also be observed in Fig.l6]in 
which it is evident that the a particles belong to a curve rather 
than a straight line. Moreover, as already stated, the observed 
c-lines are internal so that the effect cannot be ascribed to sur- 
face effects. A theoretical explanation for this behaviour and 
a detailed discussion of the energy dependence of LY and QF 
in different scintillating crystals can be found in fSOl. Essen- 
tially, the observed behaviour reflects the fact that an a particle 
scintillates less with respect to an electron because of its larger 
dE/dx which can induce saturation effects in the scintillator (the 
Birks law |21| ). The energy dependence of the LY and QF is a 
consequence of the energy dependence of the stopping power. 
Electrons, due to their low stopping power, do not suffer of sat- 
uration effects and their LY is, consequently, energy indepen- 
dent. For this reason the definition of QF„=tan(0Q)/tan(6J8y), 
given in Sec. [3] represents only an approximation. 

Using the values in the table we can evaluate the scintil- 
lation Quenching Factors of a particles {QF-LY{a)/LY(J3)y. 
QFi80n,=0.160±.006, QF2_«j; =0.192+. 006, QF2My=0.201+.006. 



Within this framework, furthermore, it is possible to evaluate 
the total amount of light k that escapes the crystal. 

Using Eq. l3] and considering that, for definition, 
LYsL/E=ySeyt, we obtain: 



LY — tan(Ofi,)ak 



(5) 



Now we have to note that the energy calibration we adopted 
for the Heat axis (as discussed in Sec. l2]i implies H=E. This 
means (see Eq. l2]| that Q;=l/(l-ky), being ky the absolute LY 
for y/yS events. Using this last relation in the last equation we 
finally get for y//? 






LY^ 



(6) 



tan(Ofi,) + LYy 

Using the values of Tab. [T] we get ky= 4.6 % 

As a final remark, we point out that the large spread in the 
light channel (6.8 % at the 2615 keV scintillation signal), as 
observed in the scatter plots, cannot be dominated by the fluc- 
tuation of the Poissonian statistics of the absorbed photons. In 
fact we have that the mean energy absorbed by the LD for a 
2615 keV interaction in CdW04 is of the order of 2.615-17.6= 
46 keV. Assuming « 3 eV/photon we get (assuming a Fano fac- 
tor=l) the variation in the light channel to be of the order of 
l/V46000/3=0.8 %. Moreover in other tests performed with 
different CdW04 crystal samples we observed the same anti- 
correlation (in terms of 9fj,) but with a much smaller spread. 
Thus the magnitude of this spread is probably dominated by 
fluctuation in the overall light collection efliciency. In partic- 
ular we believe that the grade of the surfaces of the crystal (in 
terms of diffusion of the scintillation light) plays an important 
role in this mechanism. 

6. Heat absolute scale and Heat Quenching Factor 

From the previous sections is evident that for scintillating 
bolometers the Heat/Energy scale is different for different kinds 
of interacting particles (or, better, for particles characterized by 
different LY's). In other words, with our convention for the 
Heat axis calibration, H measures the total energy of a ylfi 
particle while it doesn't for an a: the energy calibration of a- 
particles has to be dealt in a separate way. 

This "displacement" for the a lines is evident in our exper- 
imental data. In particular the '^°W-Qf line (as can be seen in 
Fig.|2|i, whose nominal Q-value is 2516 keV 122], appears (in 
the y-calibrated spectra) at 2627+2 keV. This line is rather pe- 
culiar since it definitely proves that the mis-calibration cannot 
be ascribed to a wrong extrapolation of the detector amplitude 
vs. energy calibration (the highest y-line used for the calibra- 
tion of the Heat signal is indeed the 208'pj jjjjg ^^ 2615 keV). 

The observed shift in keV, AZs - Ea - Eynj, between the 
experimentally reconstructed energy and its nominal value is 
given - for the three mentioned or's - by: 

KExmv, = 110+3 A&isy = 172±3.5 A&^^y = 180±2 (7) 

This shift can be easily recognized by looking at Fig. [3] Let 
us consider the case in which the crystal does not scintillate: 



in this case the two particles, a\ and y\ will show, as exam- 
ple, E{a\) < E{y\). If the crystal scintillates and the light 
escapes the crystal, then the two events will move left on the 
two different iso-energy lines up to reach the case in which 
H(al)=H(7l). This means that the energy shift can be writ- 
ten as h.E - //^j"" - H^'f'. By simple geometric consideration 
we get 



^Ught=Q ^^^ Lightyip ^^L_^ J^^ylP 



^ylP taniOfi,) 

So, finally, we obtain 

A£ = £ [MeV] 



J-^Yyip - LYa 



tan(Ofit) 



[keV] 



tan(6fit) 
Inserting the experimental values we get: 



(8) 



(9) 



AEn 



107+3 AEi: 



175 ±6 AEii^ 



195+7 (10) 



These last values are in very good agreement with the exper- 
imental ones in Eq.lTj confirming our model. 

To conclude this section, we need to face an important ques- 
tion about the definition of the Heat Quenching Factor (HQF). 
The conventional way to treat differences in the response of a 
detector with respect to different type of interacting particles is 
to introduce relative quenching factors. This is exactly what 
is done for the scintillation light. In the same way, following 
what is often done in literature li23l l24l |25] |26]| . it is possi- 
ble to introduce the HQF as the ratio between the Heat sig- 
nal of an a particle and that of a y/jS of the same energy, fully 
absorbed in the detector With this definition, for our detec- 
tor we can quote the HQF measured for the '^"W a decay as: 
HQFisov(/=2627/2516=1.044. It is straightforward to show that 
in our framework we can easily evaluate HQF=1+ J^'!^ ° = 

Even if the definition of the (scintillation) QF and the HQF 
is identical, due to the different physical mechanism that gener- 
ates them, they show a substantial difference. While the former 
depends, substantially, only on the "nature" of the scintillator, 
the HQF, instead, depends on the amount of light (energy) that 
escapes the crystal: it depends on the size and/or on the quahty 
of the crystal. It is also obvious that since some (minor) amount 
of light that escapes the crystal can be reflected by the reflecting 
foil and can be absorbed by the crystal itself, then the HQF will 
depend also by the setup itself. 

With the 3x3x2 cm^ CdW04 crystal measured in [9J, we 
found the position of the '^''W line at 2670 keV, giving therefore 
HQF= 1.061. 

Now one relevant question that could be addressed is: do we 
have bhnd channels in our detector? We don't have a general 
method to investigate this point, but the results presented so far 
are still valid even if we assume that a fixed fraction of the total 
energy is lost in blind channels (we used only Heat/Light ratios 
in our framework). Then we can investigate if this fraction is 
different between a y/fi and an a. In other words we could ask 
ourselves if the sum of the heat and light signals is the same for 
particle converting the same total energy in the scintillator We 
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Figure 5: Scatter plot Heat vs. Light obtained in a 8 hour measurement with 
an Am-Be neutron source. In the inset the low energy region. The energy 
spectrum extends up to 9 MeV. The 4.44 MeV "line" due to the '^C* is visible. 
The intrinsic energy spread of this peak is rather large (« 1 %) due to Doppler 
shift induced by the "in flight" decay of '-C* in the source. 
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Figure 6: Scatter plot Light vs. Heat obtained in a 93 hours ^^^^Th calibration 
measurement with a smeared a source facing the CdW04 . The liquid standard 
also contains a small amount of -'''U. The alpha curve is completely separated 
from the y//3 line. Moreover it can be noted that the a events are not belonging 
to a straight line, showing a decrease of the QF with the energy, as discussed in 
Sec.|5] 



provide therefore an alternative definition of HQF as the ratio 
ofE^''^'"=%ndE^;f'=". 

Provided that our model that compensates the energy losses 
due to scintillation is correct, using the values of Eq.[TO]we can 
evaluate the (corrected) HQF: 

HQF(i80W)=1.001+ .002 HQF(238U)= 1.000 + .002 and 
HQF(234U)=0.997 + .002. 

These HQF's are compatible with 1 meaning that there are no 
intrinsic differences in the heat signal generated by an a with 
respect to a signal generated by a yZ/Q Alternatively we can 
state that any eventually present blind channel have the same 
behaviour for the two. 



7. Background Rejection: a's and neutrons 

The detector was also exposed to a neutron source. It con- 
sists of a 185 kBq Am-Be source with a neutron production 
rate of a: 10 n/s. The neutron spectrum has its maximum at ~5 
MeV with an high energy tail reaching 10 MeV. The source was 
placed inside the shielding of the cryostat, in the same place 
where the y sources are inserted. The scatter plot is shown in 
Fig.|5j neutron direct interactions are clearly visible, especially 
in the inset. The y/yS events extends well above 2615 keV due 
to (n,y) reaction on the surrounding materials but, mainly, by 
the source itself; for each neutron produced there is 60 % prob- 
ability to produce an excited state of ^^C that emits a y of 4.44 
MeV. 

Moreover it has to be remarked (see also Sec.jSjl that "^Cd 
has a huge neutron capture cross section, with a Q-value larger 
than 9 MeV, so that "mixed events" are possible (a neutron 



'We note here that although we always discuss of o-'s it is true that we are 
always considering the system (a + nuclear recoil), where however the nuclear 
recoil transports only a small fraction of the total energy. 



scatters on the Oxygen of the crystal and then is absorbed by 
the "^Cd with subsequent de-excitation of the nucleus). From 
Fig. |5]we evaluate the neutron scintillation QF with respect to 
y/yS to be (0.14±0.03). The error is dominated by the systemat- 
ics induced by the extreme weakness of the scintillation light in 
the region Oh-300 keV, where most of the events are recorded. 

In a subsequent test, an a source was mounted close to the 
detector in order to evaluate its rejection capability. The custom 
source was build up using 2 yuL of a 0. 1 % ^■'^U liquid Standard. 
The U "droplets" were dried on a 1 cm^ Al tape and then cov- 
ered with a 6 /im Aluminized Mylar foil in order to "smear" 
the alpha energy down to the region of interest. The source was 
then faced to the crystal on the face opposed to the LD. In Fig.l6] 
we present the scatter plot obtained with the above mentioned 
source during a 93 hours ^^^Th calibration measurement. Fix- 
ing at Act the acceptance on the light signal of the 2615 keV 
y-line, we evaluate a rejection factor for a particles >35cr at the 
same energy. As far as the fast neutron interaction is concerned 
, the rejection factor will be obviously larger, but the most dan- 
gerous contribution will arise from thermal neutron absorption 
through "''Cd. But, at least in principle, thermal neutron can be 
effectively shielded while fast neutron, via /v-spallation within 
the shielding or close to it, cannot. 



8. Background measurement: internal contaminations 

In Fig. It] we present the Light vs. Heat scatter plot obtained 
in the 433 h live time background measurement. It can be eas- 
ily recognized that we observe only two transitions belonging 
to the Uranium chain, namely ^^^U and ^^^U, corresponding to 
an internal contamination of 3.1-10"'^ g/g and 5.7-10"''' g/g, re- 
spectively. The Uranium chain is broken at ^^"'U, while a con- 
tamination of ^'"Po is clearly seen. As it often happens, it is 
rather difficult to determine if this Po contamination arise from 
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anticorrelation between light and heat improves the energy res- 
olution by a factor 1.26 @ 88 keV and up to a factor 2.6 @ 
2615 keV. We developed, for the first time, a procedure that 
evaluates the Heat Quenching Factor corrected for the loss of 
the scintillation Hght. 

This CdW04, grown without any precaution in terms of ra- 
diopurity, shows extremely low trace contaminations in U and 
Th. Nonetheless a future experiment based on this compound 
should use isotopic enrichment, for two reasons: i) increase the 
mass of "^Cd; ii) decrease the fraction of "^Cd in order to 
avoid y6 decay and, more important, decrease (n,7) reaction due 
to'i^Cd. 

As a final remark we point out that such a crystal, depleted 
in "^Cd, would be an extremely interesting compound not only 
for DBD experiments but also for Dark Matter searches. 



Figure 7: Scatter plot Light vs. Heat obtained in a 433 liours live time back- 
ground measurement. Part of the ^'"Po "bump" is characterized by smaller 
light emission. This could arise from a worse light collection efficiency, espe- 
cially from the lateral part of the crystal. The energy scale of the a is corrected 
according with what exposed in Sec.l6] 



its parent, Pb. Moreover the absence of a clear peak indi- 
cates that the contamination is partially penetrating (few fim) 
the surface of the crystal (or its surroundings). No events can 
be ascribed to the Thorium chain, giving a limit of Q-IO"'-' g/g 
(95% CL) for ^^^Th, one of the most "dangerous" contaminant 
for DBD searches. But the main (unexpected) feature result- 
ing from Fig. It] is the y/yS background above 2615 keV. These 
1 1 events cannot arise from external y's due to ^os-pj qj. 2i4gj 
since the observed background lines (583 keV, 2615 keV and 
1764 keV) are too weak to allow such contribution at high en- 
ergy. The same holds for internal contaminations belonging to 
Th and U chain. Unrecognized pile-up of the "''Cd-yS decay 
(Q= 320 keV) with a 2615 keV y cannot contribute at this level 
(we expect ~ 2 ev/year). Internal contaminations due to rare 
high-energy y6 emitters like '"^Rh (Q=3541 keV) are extremely 
difficult to evaluate. But, in any case, we observe j//3 events 
up to 4.5 MeV, extremely large to be associated with standard 
"known"/? emitters. On the other hand we previously tested dif- 
ferent CdW04 scintillating crystals [281 . obtaining, with larger 
statistics, no events above 2615 keV. The main difference with 
respect to fSSl is that in the present work the environmental 
neutron shielding (consisting of 7 cm of polyethylene and 1 cm 
of CB4) was removed from the top of the cryostat, leaving 
~0.6 m^ opening, corresponding to a; 8 % of the total neu- 
tron shielding. We think that the thermal neutron absorption 
by "^Cd is the most probable mechanism in order to explain 
the anomalous background above 2615 keV. 



9. Conclusions 

For the first time a "large" scintillating bolometer to be used 
for future DBD searches was fully characterized in details in 
terms of energy resolution, internal contaminations and parti- 
cle identification capabilities. It was shown how the use of the 
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